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I. INTRODUCTION 


Tt\« thaory of th« formaclon of t;h8 solar chromosphara and corona has as 
Ics main challanslng goal Co spacify Cha natura of »na hasting mach&nitm. 
Since cha discovary of tha U I0‘K) hot solar corona about forty years 
”.,w "-wory cure".-! UiJAsini haj vc'4l**id *.*i.»Ld Lntar^'sit:. 

Chaorias ware assantially based on the dissipation of acoustic waves or 
weak shock waves generated by the photospheric grastulaCicn (Biermann, 19465 
1948; Schwarsschild, 1948; Schatssman, 1949). They looked very much like a 
description of th« extension of the photospheric model atmosphere by in-* 
troducing in addition to the radiative energy flux and Che opacity a 
mechanical energy flux and a wave absorption coefficiant. The atmosphere 
was assumed to be largely spherically symmetric and all the model 
calculations were one dimensional and steady, specifying boundary con- 
ditions at appropriate levels in the atmosphere. 


It was soon realized chat the corona overlying an active region needed to 
be heated one to two orders of magnitude more efficiently than the so-call- 
ed quiet corona and because of the active region's magnetic nature it be- 
came evident Chat magnetic fields must play an important role in the 
energetics. 


Av.Jt.v.r important discoverj* was the solar wind which is a direct con- 
0 -* the e:dscence of a hot corona and thus of coronal heating. 
Spacebome observations of the solar wind made it clear that actual 
irreversible heating as well as momentum addition must be present up to 
large distances from the sun, far beyond the level of the maximum 
coronal temperature. 


It is unlikely chat weak shock waves can reach that far. They are probably 
exhausted already in the inner corona. Here the magnetic field must play a 
dominant role in selecting Che required mode of energy transfer, p: ».suraab- 
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ly in tho form of AlfvSn wav«$i, 

OacASlonally ch« corona ia haacad Co vary high CamparaCurat (lO’^ “ 10** K) 
4 « ; *141*4. Tt is tauar.nUv '•-oeapcad that Joula haating and 
magnstic raconnoctioni convarsion of s corad masnacic anargy Into Uanc and 
high anargacic parcialaa, is the fundanancal flara procasa, and wava 
dissipation plays a nagligibla role in this case. 

t^a are chus facing cwo basically dlffaranc mechanisms of coronal haating; 

a. Heating by acoustic processes in the ''non magnetic" parts of the atmo- 
sphere [ S ■ ( 8irp/B^) Z I } of which the shock wave theory is 

an example worked out in great deal. 

b. Heating by electrodynamic processes in the magnetic regions of the 
corona. ( S « 1 ) either by 

0 magnetohydrodynamic waves or 

2} current heating in the regions with large electric current 

» 

densities (flare type heating) . 

The pr«S"?nt state of the art in rhase two classes of theories and the 
typical problems one encounters in each of them will be discussed in 
sections 3, 4. 

,u «; 441144 Che shock wave theory was rather successful in that it linked 
the observed phetospheric motions to the observed and predicted chromo- 
spheric and coronal energy losses known at chat time. However, the amount 
of mechanical energy generated in the sub-photospheric layers is known 
with an accuracy less than an order of magnitude and so are the radiative 
losses of Che chromosphere and the corona. 

For an extensive discussion of models based on shock wave theory the 
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rsadav is cafanred to n great number o£ $;aviev) papers (Kuperus, I9ft9; 
Bray and Loughhead, 1974; Kaplan at al,, 1974; Athay, 1975; Kuperus and 
’.’icUHria nn-l ’'’■>■'^ 1 . ''^77! yiowt and Pineaii de» PorSts, 
1977; Ulmschneider, 1967, 1970, 1971). 


At present tha theory of the fiv>naation of the outer solar atmosphere 
should take into account the following reasonably well established 

facts. 

I, The corona consists of essentially two types of radiation structures 
shich could be associated with or coincident with magnetic 
structures, (Vaiana and Rosner, 1978), 

a) OPEN’ STRUCTURES, the so called coronal holes, which are thought 
to be the seats of high speed solar wind streams and have a very 
low emission measure in soft X-rays (Eirker, 1977). 

b) CLOSED STRUCTURES (loops) that radiate more energy than the 
surrounding corona. 

The loops are strongly concentrated in the active regions but 
also interconnecting loops are present between several centres 
of activity. 

.. * ..^jp’.Yw’cic network and the plage regions are the seats of 

strong magnetic fields (~ 2000 Gauss). The field appears to be 
concentrated in flux tubes in the lower parts of the atmosphere. 

They appear to have the same structure underneath the coronal holes 
as elsewhere (Stenflo, 1978; Zwaan, 197S; Spruit, 1980). 

3. The existence of a thin transition region between the chromosphere 
and the corona with steep temperature gradients and large "turbulent” 
velocities confirmed by the XUV observations (Withbroe and Noyes, 
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1977; Bruecknsr, 1980). 

4. A highly filamencary chromoiphara, with upward motions in tha 

rad dovnflows near tha network boundaries (Brueckner, 

1980). 

5. A photospheric/chromcspheric velocity field with turbulent properties 
which apparently are evanescent acoustic waves chat have a spectrum 
peaked around a period of 300 sec. (Deubner, 1977; Stein and I.Qih<* 
bacher, 1974). 

6. The detection of a strongly varying multicomponent solar wind plasma 
at Che Earth's orbit with large differences between Che electron and 
proton temperature , a very distinct magnetic sector boundary struct-* 
ure and some evidence of the presence of AlfvSn waves (Belcher and 
Davies, 1971). 

In the next section the requirements of chromospheric and coronal 
heating will be discussed in the context of the fundamental constraints 
one encounters in modelling the outer solar atmosphere. 
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11, ENERGY TRANSFER IN THE CHROMOSPHERE AND CORONA 


lli l, Funcksrgntal SqtuxHons and Brntp-jy B^Uxma 


ModalUnx, chs oucor *oXat atmosphere involves solving the equations of 
continuity, mymantuia and energy in a given geometry usually specified 
by the nuignetic field ^ and the acceleration of gravity . The 
continuity equation is - 7 . p v , where p is the mass density 

and ^ the velocity. The equation of motion is 

s1| 4 - J 

o C ^ 

where p ■ RpT/u is the gas pressure, T is temperature, R is the 
gas constant, y the mean molecular weight, g - -7i{) , where is 
the gravitational potential, ^ the magnetic field vector, 
j is the current density and T is the viscous 
stress tensor (Landau and Lifshits 1959). 

The equation for the thermodynamic energy U ■ i P + pe + pi) , 
where e - p / (Y“I) is the internal energy, is 



where the total thermodynamic energy flux F is given by 






(3) 


■RAD 


The terms on the right hand side are respectively the convective 
energy flux F , the mechanical energy flux F , the viscous 
energy flux , the conductive energy flux » ®nd the 

radiative energy flux Xrad ‘ ^ coefficient of the thermal 

conduction, which will be discussed in the next subsection. 

The term j.E is the electrodynaraic energy input representing Che 
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coupling bitcwccn cho Chtnsodynatnic tyititm nnd cht *l»ccrodynnmic 

«y«eow. E * nj ”■ Cv X B) g”* it the total •Icctclc field und n 

the aluccri.'al rasiafeivity. Inataad of writing the anargy aquation for 

cue cuv.raucj .in»aie anacg) one couW write thw equation for the total 

anargy including the wagnatic energy* Using Poynting's thaorem 

the energy equation can be written in conservative fora 

3U /3t ■ -V.? , by adding the Poyncing flux (E x B)c / At to the 
total flux. We prefer to write the anargy aquation in the fora of 

aquation (2) clearly separating the alectrodynaaic coupling tara 

J * 


Viscous coupling in the outer layers has not been considered 
since small wavelengths are required for the viscosity to become 
important. The viscous flux will be neglected further. In a non 
magnetic atmosphere in steady state the energy equation (2) is thus 

0 . 

The energy requirements for the heating of trie chromosphere and 


” ■ ^Gorvv .imech «~cand «**rad^ 


' . . .r ' Hued on the strongly model dependent estimates of the 

Cvergy looses in radiation, conduction and convection. Extensive 


discussions of this problem are found in Athay (1976) and Withbroe 
and Noyes (i977). We here confine ourselves with summarizing the total 
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•gtimactd tnorgy loisti of th« chrono«phtr« and corona in eha quite 
gun I cht coronal holt end cht tceivt rtgion at hat bttn given by With- 
brot and Hoytt (1977). 

Thtit tntrgy fluxtg havt eo bt rcpltnightd by acougeical and/or 
.■iiuscodyn i®ic heacir-i* A fluK Ir. batwaetv 10* ~ !0’ ergs* / cm* huc 

is ntctggary dtptnding on Cht dtgrtt of activiey and cht nagnttic 
fcopology in cht corona. Tht tntrgy dtmands art Cht gmallttC for Cht 
quite gun alchough Cht rtdiacivt losgtg art an ordtr of Mgnieudt 
largtr *:han in coronal holts. Tht largt tntrgy flux rtquirtd for a 
coronal holt tggtncitlly rtsulcs from cht largt convtccivt lotsts ouc 
of Cht presumably optn magntcic scrueeurts. Boch cht quite sun and 
accivt rtgion consise of closed scrueeurts and chus Che losses art 
primarily radiacivt and conducCivt. 

II. 2. hEoTOsaopiQ Ef feats and Tvanapovt Fapametavs 

Wt ctnnoc over emphasize the care one needs Co exercise when one con- 
siders Cht tntrgy and momtncum balance problem in a magnetized plasma. 
This arises for four reasons: (i) Cht magnetic field can result in 

highly anisotropic transport coefficients; (ii) some of the ED heating 
mechanisms cause localized heating which can be isolated by the 
ra^inetic field surrounding the heating locus; (iii) some of Che ED 

ruBchanisns can lead to substantial modifications of the trans- 
port coefficients themselves - at least locally - and these modifica- 
tions must be accounted for self-cons is tently when treating Che energy 
and momentum balance question; and (iv) some ED mechanisms are long 
wavelength mechanisms (a mechanism that dissipates energy over a sub- 
stantial fraction of the magnetic configuration parallel to while 
some ED mechanisms are short wavelength mechanisms (a mechanism Chat 
is highly localized parallel to (Spicer and Brown 1980) . Such 

( original vagr is 

Of QUALITY 



behaviour on bhu parb ol tht htating iMchani»iat can iigniCictnbly alccr cha 
•ncrgy balanci aeguninbi and caauleing lualing lava (Spictr ab al 1980). 

•{ bho haacing machanlam opacabaa along bhc magnabic field over a langbh 
L|j ^ I vhich i i bha charaebariibic langbh acroaa bha ntagnabic field, 
,»aau .K % uia i»c the divcrg<‘nc » of tiid jwai. fluj; paralUl end parpandlcul- 
at bo bha magnabic field can be of bha order one. Such a tibuabion appliea 
bo bha bearing node proeaat (Spicer 1980) (IV). in the caie of anomaloui 
Joula hearing bhia rabio may be much amallar bhan one. 

The above argumanba should be conbraabed wibh che classical braabmanb of 
bhertnal conducbion (Chapman 1954; Braginskii 1965; Kopp 1968). Thera ib is 
shown bhab bha tieab flux parallel bo B is given by 

-vcond II " " ^0 ^11 ^ * where T*^® • 6 x 10”’ T*/* erg cm"‘ 

sac”‘ K”^ . This flux is primarily direebed along JB , unless ' • 

This uas firsb poinbad out and applied bo magntbic sbrucburaa in the 
chromosphere and corona by Kopp and Kuperus (1968) and labar reexamined by 
Gabriel (1976). 

The classical apptoxiaabion neglecbs che cr>llaccive flucCuacions. 
However, when a non'-bhermal faabure imposes ibsclf on a plasma che 
coUecci’.’e fluccuacions can be driven co higher amplitudes, These collect- 
ive fluctuations can then produce large stochastic particle deflections 
• hus enhanca che effective collision frequency thus leading to 
transport parameters. The non-themal feature that is of 
present interest is the drift velocity of the current (IV’) . which must 
reach a magnitude comparable to the phase velocity of a normal mode of 
the plasma. 

Another effect that is microscopic in origin but also manifests it- 
self macroscopic«illy is the divergence of the radiation flux , 

which can be hoth a .sink and a source term. This point is generally 
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ovtrlook«d in «ntr$y balanca a«gum«nfcf Withbroo nnd Noyti 1977) 

•van chough tffloreicanc haaCing of tha chromoaphatra by * I kaV goffe 
S-rayg aniccad in the corona by fXacag is known to significanely altar 
tha anacgy balanct arguiBunti in solar flares (Uanoux and Nakagawa )977{ 
Uanoux and Ruse !9d0). Similar affucca could also occur during ((alascont 
periods and could radically alcar cha prasanc anargy balance arguments 
(a.g., Uithbroa and Koyas t977). 

Houaver, since the corona is Co a large degree optically thin one may 
use cha radiaeiva loss funccion ■ p'* f(T) , where f(f) has 

bean cclculacad by Cox and Tucker (1969) and lacar refined by McWhirtar 
at al 0975). 

An addicional microscopic difficulcy associatad wich soma ED 
mechanisms is cha problem of how energy ralaasad in a localised volume 
dV is rtdiscribuced chroughouc Che cocal volume V of Che magnecic 
configuracion such chat dV/V « I . Physically che problem, is to 
decermine which cransporc mechanism is best suited to Cransporc energy 
perpendicular co B from the locus of energy release inco che 
remaining volume V . Is ic diffusion or convection? Diffusion process- 
es are rather unlikely due to the large timescales compared to the 
characteriscic cooling timescales. This letter point is true even with 
’.••• ernnssort. However, convective transport which leads to much 
ihorcsr cimescaloi can accomplish this by a variety of means usually 
depending on the ED heating mechanism. For example, tearing modes (IV) 
could lead to magnetic braiding as has been discussed in connection 
with flares (Spicer 1976, 1977) and thereby lead to v». substantial 
increase in the heated volume over the actual energy release volume 
dV . On che other hand local energy release roust lead to waves being 
generated at the site of energy release. These waves can then propagate 
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and damp throughout the volume V thus heating the plasma as they 
proceed. This approach has been cKamlned by Spicer et al (1980). 

11.3. Global Slnopgy Balance of'td l4aa)*08aopio Effects of the Magnetic Field 

11.3.1. GL0B.\L ENERGY BiVLANCE AND NON STATIONARY HEATING The formation 
of coronal structures requires a detailed knowledge of the local heating 
throughout the whole structure. Nonetheless several global calculations 
have been presented balancing the energy input against the total losses 
in order to find the dependence of the coronal temperature on the energy 
input. 

Rosner et al (1978a) derived a scaling law for coronal loops: 

T - (pL)^/^ , where L is the loop length and p the pressure in the 
loop. Hearn (1975, 1979) suggested that at a given base pressure and a 
given temperature the corona will adjust in such a way that it picks out 
the solution with a minimum energy flux, thus reducing the number of 
parameters. 

However, Endler et al (1979) demonstrate that it is the details of 
the heating mechanism expressed for instance in the dissipation length 
chat determine Che ultimate coronal configuration. For instance the 
height where Che maximum temperature is reached strongly depends on the 
• ; (Vlmschncider 1971; Mchlairter et al 1975). The basic 

problem with the mininium flux assumption is that it seems to be over- 
determined. Once Pg has been selected the radiative losses are 
essentially decermined and so is the energy input. Changing tho energy 
input will result in a shift of the transition layer to a pressure 
such that energy balance is maintained (Martens 1980; Van Tend 1980). 

How this occurs can only be found after solving the time dependent equations 


ORIGINAL PAGE IS 
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Ev«n in th« clmpU cnit of plan* gaoroefry iuoh haa not ync bean done In 
any *aci*faccoti» way* Kuperue and Athay (1967) at|uad tnat a jeiry aceep 
iowei* ttranaicion layav ia vairy hard to heap in hydcotcacic eqniUbrinm. 
rh# ar$ufflenc is thac cHa himh coronal camptrature ia a natural con- 
aequanca of the inability of the upper chromosphure to radiate Che 
disiipated energy in a low temperature configuration. IE it were not 
for the thermal conduction two eeparated but adjacent atmospheric 
layers could exist both in hydrostatic equilibrium aasumlng pressure 
continuity. However, the coefficient of heat conduction very abruptly 
decreases going inwards from the transition layer to the chromosphere. 
Tlterefore, if thermal conduction is the essential energetic link 
between the corona and the upper chromosphere, the upper chromosphere, 
which ia supposed to radiate at its maximum efficiency ia likely to be 
excessively heated. As a consequence the upper layers may be set into 
motion* Kuperus and Athay thought that this overhaating in the upper 
chromosphere would be the origin of the spicules. It was demonstrated 
by Kopp and Kuperus 0968) that overheating aggravates in converging 
magnetic fields due to the magnetic chaneiling. Converging flux also 
increases the likelyness of heat flux instabiUty (Spicer 1979) , This 
would then argue in favour of seeing the spicules concentrated on the 
'• ndu'ies whei-e one finds the magnetic knots. 

li, i.! tempting to assign the observed large non thermal motions in 
the transition layer (Bruacknar 1980) to the effect of dynamic adjust- 
raent of the transition layer to time dependant coronal heating. 

On the other hand Pneuman and Kopp (1980) showed that the observed 
downflows in the transition region carry an anthalpyflux larger chan 
Che conductive flux so that in those regions the models based on con- 
ductive energy supply need to be reevaluated. 
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Static modalling ia Ukaly to b« of Uttla vjiluo and miaht; avan lend 
to arronnoua conctuaioas . 

Uaaaattly tha magnitudia of tha tamparatuta gradiant in tht tranaition 
cone has b««n questioned on thaotatical grounds (Spicer 1979) and on 
arapicicai grounds (HicoUa at al 1979). 

Using empirical results Spicer (1979) has argued Chat because the 
mean free path of Che electrons in tha transition zone is comparable to 
tha empirically determined characteristic cemperaCura gradient scale 
length there the heat flux will be marginally unstable electrostatically 
if the empirical models ate correct. He showed that tha heat flux would 
be reduced by - 30 over its classical value. Spicer further argued 
that this reduced heat flux should appear as convective turbulence in 
the transition zone consistent with the previous arguments of Kuperus 
and Athay (1967). 

On the ocher hand Nicolas et al (1979) using high resolution d.ata 
( ~ 1" ) from the HRTS instrument (Bartoe and Brueckner 1975) have 
brought into question tha assumption of taking the filling factor equal 
to I for slit instruments. They conclude by suggesting chat most of 
Che transition zone emission may origi. ace from structures whose 
filling factor would be less than one, hence reducing Che temperature 
• «• *.<nficCad. 

II. 3. 2. HACROSCOPIC EFFECTS OF THE IIAGNETIC FIELD As already discuss- 
ed Che magnetic field ^ plays a crucial role in the coronal hearing 
process. However, it is not yet clear whether this role Is passive, 
chat is B simply acts as a guide for the transport of energy and 
momentum, or whether this role is active, that is B itself is in- 
timately involved in the heating process. This delineation of roles for 
B naturally separates the two classes of coronal heating mechanisms : 
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acoustic (III) and alectrodynamlc (IV). Howaver, such a pracis® d®*- 
lintntion of toUs for Q is quit® artificial sine® for h. not to 
play a rola in tha macUanical class of coronal haating machanisms ra- 
quiras tha acoustic wavas involved to propagate exactly parallel to B , 
since oblique propagation will lend to tnagnociu field perturbations 
iS^ that land to currants and ion oscillations that damp by Joule and 
viscous dissipation. Since tha generation of 5B. represents work done 
by the acoustic waves oblique propagation naturally reduces the damp- 
ing length generally assumed (cf. Ill) of the acoustic waves. Neverthe** 
less, we will assume this delineation to be valid in the zeroth 
approximation. Under these circuq>8tances, B plays an active role only 
in ED mechanisms and this automatically imposes strong constraints on 
both the topology of B and its magnitude. 

As we will discuss in IV three ED mechanisms (surface waves, re- 
connection and anomalous Joule heating) require relatively large 
gradients ir B and lead to heating in narrow layers. Since large 
gradients in B implies large current densities we are faced with not 
only understanding the ED heating mechanisms but also understanding 
the global stability of magnetic configurations in which some of these 
machanisms are supposed to operate. 


r-..cnt observational limitations, little can be said of the 
actual internal structure of the magnetic configurations observed. 
However, there are few theoretical points chat can be made. 

The first point is chat steady state heating by definition is im- 
possible in force free fields (J II B) , because the dissipation of 
magnetic free energy must be balanced by a divergence of a Poynting 
flux, which by definition is aero for a force free field. Hence, 
coronal heating by an ED mechanism which operates in a force free 
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field, will b« cioe dependent because it must occur in a time at least 
shorter than the radiation time scale. 

The second point is that any ED coronal heating mechanism that 

dissipates magnetic free energy must heat the plasma in a time t^ 

shorter than the characteristic plasma cooling time T . This Implies 

that the growth time t^ of the mechanism must be shorter (generally 

t < t ) , which for the anomalous Joule heating and reconnection 
g q 

mechanisms requires strong magnetic fields and/or strong magnetic field 
gradients (cf. IV). However, strong fields and field gradients are 
difficult to maintain in the open corona without external stresses due 
to the tendency of the field to expand so as to relax the gradients and 
weaken the field. Hence, to neat the corona either in steady state or 
quasi-steady state by an ED process such as anomalous Joule dissipation 
it is necessary for the divergence of the Poynting flux to resupply 
energy in a time shorter than a plasma cooling time (by conduction, 
radiation or expansion) and also with a magnitude sufficient to 
compensate both for energy dissipated by the ED mechanism and the 
energy lost to expansion. Such constraints impose further constraints 
on the current generators which are necessary for ED mechanisms to 
function. 

*.•. principle the whole current system can be defined by a knowledge 
of the velocity field, magnetic field and the pressure gradients. This 
can be seen by deriving ^ from the momentum equation. Neglecting 
gravity one obtains: 





Cf df<y -5 
, £. a. X 3 

3 ^ cLt ^ 


9 


C 


( 4 ) 
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From V.j ■ 0 one findi ji. after integrating along a field lines 



Thus a solution to the appropriate equations defining ^ , B and p 
together with the observable boundary conditions on ^ ^ P 

the photosphere, do specify ^ ^ and p in the solar atmosphere and 
thus j|| and . 


(5) 
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III. HEATING BY ACOUSTIC PROCESSES 

The ehtory chat haa baen popular lor mora than 20 yaars la that tha 
chroDosphara and tha corona ara haatad by pariodic ahucK wavaa ganaratad 
by tha turbulant motions in tha uppar layars ol tha convection zone. Tha 
theory is based on the pioneering work ol Biemann (1946), Schwarzschild 
0948) and Schatzman (1949) and later extended by Da Jagar and Kuparua 
(1961) • Ostarbrock (1961), Uchida (1963), Bird (1964), Ulmschnaidar 
(1967), Kopp (1968) and Kuparus (1965, 1969). 

Tha basic idea behind this theory is that due to turbulent pressure 
lluctuations a liald ol small amplituda acoustic waves is generated. 

Due to tha axpotientlil dacraasa o£ tha density the velocity amplitude 
o| the waves as they travel outward incruasas rapidly according to 

As tha velocity increases non-linear ellacts of the propagation causes 
the sinusoidal wavs prolila to steepen into a saw tooth like sequence 
ol shocks. In these shocks the dissipation rate is much larger than in 
the linear acoustic waves. Moreover the heating rate is thus simply 
.ecarsin^aj by the shock wave Mach number, , and the period P . In 
this section we will first discuss the generation ol acoustic waves 
and Chen the propagation and dissipation ol shock waves. 

1. Gemmation of AoouBtia Wavss 

Acoustic waves are generated in the upper convective layers of the Sun 
by compressional disturbances associated with the turbulent motion. In 
order to calculate the total amount of acoustic energy leaving the 
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convection 7.one iolar physicists have profited from the very elegant 
theory of the aerodynamic generation of acoustic noise by turbulence 
developsd by Ughthill (1952, 1954, 1960, 1962), 


Tf 0, is the density fluctuation it follows from the equation of 
motion and the continuity equation that 




(7) 


where i,j ■ 1,2,3 indicate the three coordinates x,y,z, v is the 
sound velocity, are the forces acting on a fluid parcel and S. . 

* J 

are the turbulent stresses. F. ■ p, g 5; and S.. ■ p v. v. . 

l » is 11 , I j 


It has been shown by Lighthill that there are two equivalent ways of 
considering the turbulent field. In the first place one may consider 
the turbulent density p^ at every point as a consequence of the 
internally acting forces, F{ and stresses, S,. . Secondly one may 

* Ij 

consider the turbulent region a whole subjected to an externally 
applied system of fluctuating forces and stresses. In both cases 
equation 7 describes the generation of sound. The sound producing 
sources are at the right hand side. According to the way the forces 
and stresses are differentiated we distinguish between dipole 
r (she first term) and quadrupole radiation (the second term). 

The intensity of the sound waves is given by I(x) ■ p^"' (x,t) > , 

where the brackets denote an averaging over time. The total power out- 
put due to dipole radiation is - 6 c and due to quadrupole 
radiation ■ a c where e ■ p < v^ L“‘ is the turbulent 

dissipation, L is the characteristic turbulent length (e.g. the 
mixing length) scale, M is the Mach number of the turbulent motions, 

M ■ < v^ v”' , and a and S are constants depending on the 

s 
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eurbuUnt correlations a ■ 0(10) , 0 ■ 0(0.1) . Hence Pq/Pj ■ 

• a S'* M® >1 for M > 0.1 (Unno and KaCo, 1963j Unno, 1964} 
Kuperus, 1965). 

Magnetic fields in the turbulent layers may add to the dipole emission 
(Osterbrock, !96l) as well as to the qaadrupole emission (Kuiscudi 
1955). An essential difference between dipole and quadrupole radiation 
is that dipole radiation is strongly peaked in the direction of the 
forces while the quadrupole radiation is isotropic. According to Kuls- 
rud (1955) the amount of noise generated in the magnetic regions 
(8 ** 1) may be an order of magnitude larger than in the non magnetic 
regions (0 » 1)*. The total flux of energy (F||^) can be found by 
integrating the acoustic power over the depth. Fg, ~ 10^ - lO' 
ergs/cm‘ sec for quadrupole radiation. The uncertainty is essentially 
caused by the uncertainty in the turbulent velocities (Kuperus, 1965). 
Remember that P ~ v* ! , an increase of lOZ in v results in more 
than a factor 2 in P . 

It has been shown by Kuperus (1972) that this emission can be con- 
siderably enhanced in turbulent convection. The quadrupole emission of 
a turbulent eddy moving with an average convective velocity v^ is 
increased by a factor ( I - v^/v, J in the direction of v^ 

oy the fact that a point overlying the noise generating layer 
"sees" a relatively larger acoustic eddy when it moves towards it than 
when it moves away from it. 

Moreover, due to the Doppler effect and Che existence of a gravitation- 
al cut-off frequency u ■ Ys/2v for sound waves the transmitted part 
of the acoustic spectrum is strongly modulated by v^ . These two 
effects are taken into account in the above mentioned amplification 
factor. As a conclusion we may state chat it appears likely that a 


iufCici*neI)r Urg« llux of acouitio w«v«s is gsntraftd in th« lurbulanc 
layers of tha upper convection zone and this flux is strongly varying 
over Che surface of the sun depending on the convective state of motion 
as veil as the magnecic structure (iupergranulation and network 

boiindariea^ . 

Moreover it is expected that the part of the mechanical energy flux in 
high frequency sound waves is strongly time dependent since it can be 
modulated by low frequency oscillations in the upper convection zone 
and photosphere. 

Tttis spatial and temporal behaviour of the acoustic energy flux must 
be reflected in the structure and development of the overlying layers. 

The existence of a gravitational cut-off frequency results in a trans- 
mitted and a reflected past of the turbulently generated acoustic waves, 

A second cut-off frequency exists the socalled Brunt-Vaisala frequency 
CDj • (y-I) g Vg“* which is slightly smaller than . 

Waves with u < oi, can propagate as internal gravity waves with small 
compressive properties in the convectively stable parts of the atmo- 
sphere. Inside the convection zone only pure acoustic waves can be ex- 
C'oore and Spiegel, 1964; Moore, 1967). However, the gravity 
vavea can be excited by overshooting from the convection zone. The 
spectrum of these waves has been calculated by Stein (1968). Both the 
acoustic and the gravity modes are refracted in the upward direction 
because the temperature first decreases. After the temperature minimum 
is reached, the waves have the tendency to refract away from the radial 
direction so that a possibility exists for trapped waves in the chromo- 
sphere. Uchida (1965, 1967) demonstrated that the transmissivity of 
gravity waves is several orders of magnitude smaller chan for acoustic 
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vavts. Indtcd th« ob»«rv«d oscillatory character of the photosphere 
confirms the existence of evanescent wave modes peaked around a period 
of 300 sec. with strong acoustic properties (Deubner, 1975 j Stein and 
Leihbacher, J976; Ulrich, I977j Cough, 1978), while little energy is 
contained in the gravity modes. Moreover It has been suggested by 
Souffrin (1966) and Schatzoan and Souffrin (1967) that the low frequency 
gravity modes are strongly damped by radiation. Xt is thus unlikely that 
these waves play a significant role in directly heating the corona as 
has been suggested by Whitaker (1964), although they may be present in 
the lower chromosphere. It should be mentioned here that due to the 
rotational character of these waves an interesting possibility exists 
for coupling with Alfvdn waves (Lighthill, 1967). In that case the 
magnetic field may pick up part of the gravitational wave energy and 
transport it with much less losses through the chromosphere and the 
transition region. The propagation of magnetohydrodynamic waves will be 
discussed in the context of the electrodynamic processes. 

2. Propagation and Dissipation of Shook Wavas 

As the acoustic waves travel upwards their amplitude increases due to 
the decreasing density. Therefore non linear effects occur deforming the 
The velocity of a point in the wave profile depends on Che 

particle velocity v . 

Because of the difference of the velocities between the crest and the 
trough of a sinusoidal wave, the distance traversed by a wave of wave- 
length X after which the crest and the trough coalesce is 
d ■ VgX / 2v (y+I) . It can be shown that shock waves are formed after 
three to four scale heights (Kuperus, 1965; 1969). 

During the propagation of a shock wave through the atmosphere some of 
the shock wave energy is transferred into thermal motions due to 
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viscous proctssas. This shock wav# dissipation is a ditact consaquanca 
of tha anttopy jump through a shock transition. This antropy jump is 
uniqualy doterroinad by tha shock wava Mach nuroh'^r or by the shock 
strength S - dp/pj where dp is tha prassura jump and Pj the 

t^reosara in tha prasviock medium. 

In a quas is Cationary atmosphara parmaatad by a saquanca of shock vavas 
separated by the accompanying system of expansion waves the energy rise 
after a shock transition is balanced by non adiabatic processes such as 
radiation and heat conduction during the expansion. 


It then becomes relatively limple to calculate tha heating during one 
cycle as has been outlined by Schatzman (1949). The amount of heating 
during such a quasi stationary cyclic process 




( 8 ) 


is immediately connected to the net increase in enthalpy, h , after tha 
medium has been expanded to its preshock conditions. The net heating 
race is then - p ^h / B . For weak shock waves Che heating 

function given by 







(9) 


In general the enthalpy jump is a complicated function of given by 
the Rankine Hugoniot relations. The propagation of a sequence of shocks 
in a stratified medium has only been treated in an approximate way. 
(Kuperus 1965} Ulmschneider 1967, J970, 1971, 1974; Kopp 1968; Flower 
and Pineau des Forets 1977). 


If one assumes zero net’-displacement, thts energy in a shock wave D 



2a 


can bi wricctn ai 

C )*1 ( 10 ) 

where brackaei denoce cha juop chrough eha ahock front and 




(in 


is tha socallad Cima'-anargy incagral (Oonezi ac al., 1977) » Drinklay 
and Kirkwood (1947) calculata cha shock propagation by assuming t^ 
proportional to tha charactaristic tima for tha vaka 


•bt 

In this casa axact diffarantial aquations can ba obtainad for dt^/dx 
and for dO/dx . An immadiata eonsaquanca of tha Brinkley Kirkwood 
thaory is that dt./dx > 0 . which isaans that tha waka of a shock in- 
craasas as tha shock propagatas. It thus appears that if shock waves 
are generated periodically (e.g. by staapaning of a sinusoidal sound 
wavs with period P), a subsequent shock catches up tha wake of the 
preceeding shock, so that t^ ** F and therefore shocks cannot be 
considered a« separated shocks. 


Actually each shock evolves in the wake of the preceeding shock. In 
... it is customary to assume that t^ i* constant so that one 
;;an derive an implicit relation for the shock Mach number 

^ (13) 

where and p^ are the density and the pressure of the undisturbed 
atmosphere through which the shock train travels. 


Another way to treat shock wave propagation is to replace the atmo- 
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sph«r« by a griac numbur of discr«C« laytiri and ealculaca cha brans** 
nisiion through aaeh contact dltcontinuity. This is the Chisntll 
sMtbod (Chisncll. 1955) that is only applicable if the reflected 
shocks are neglected. 

TUa titlru matUod tuac is fruitfully applied is the soaalled t^iitham 
method C!4iithaa, 1957) which is a modification of the method of 
characteristics and essentially means the substitution of the shock 
relations into the equations of motion in characteristic form applied 
immediately behind Che shocks. It is essential for Whicham's method 
to apply the characteristic relation along the slightly different 
shock direction. Intuitively it comes close to neglecting the 
negative characteristics. The accuracy depends on the degree of 
curvature of the positive characteristics. Extensive discussions 
of tifhicham's method as well as exact expressions for the shock 
strength gradient are given by Stefanik (1969) and Kopp (1968), who 
applied this method to the heating of the solar corona. The reader 
is also referred to Gonezi et a|. (1977) who question Che validity 
of this approximatiou particularly in the case of periodically 
generated shock vifaves for which t^ » P . Gonezi et al. (1977) give 
a self consistent formulation < a shock heated atmosphere, tor a 
niJ***k :-itrength S the maximum flux consistent with S is 
rca.aci wii«ia the shocks overlap. On the other hand a given flux of 
shock waves, F , determines a minimum value for S . Applying this to 
the upper chromosphere results in strong shock waves (S > 10) . 

Their relation betw^i^n F and S is 



where p and T are the average pressure and temperature. It follows 
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£rom (9) that in tha lower transition region where the temperature scale 
height is much less than the dissipation length for the shock waves so 
that in a first approximation F constant the shock strength rapidly 
decreases with height until S I . This is a direct consequence of the 
reflections of the waves against the temperature gradient. The atmo- 
sphere is thus roug' , divided into three layers: 

1. the chromcsphere where the acoustic waves steepen into shock waves; 

2. the upper chromosphere and the lower transition region which are 
permeated by strong shock waves. 

3. the upper transition region and lower corona which could be heated 
by weak shock waves though electrodynamic mechanisms are likely more 
relevant as will be discussed in IV. 

It should be noted that the upper chromosphere and the lower transition 
region are Just the sites of the most violent motions. (See II) 

Large amplitudes are confined to a narrow dynamical boundary layer. 

Tlie acoustic theory seems to be applicable to the lower parts of the 
outer atmosphere, e.g. chromosphere and transition region outside the 
magnetic concentrations. 
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IV. HEATING BY ELECTROOYNAMIC PRC CESSES 

IV. I . Opevoiew 

O” n*-’iCTRODVJ!A:-JIC EMEROY EUctcodynaaic processes 

can play a critically active role in Che upper chromosphere » transition 
layer and corona where supposedly 3 < I as well as in the magnetic 
concentrations throughout the whole atmosphere* (J. Hollweg, 1979). 

The generation of electrodynamic energy in the solar atmosphere 
arises when plasma is forced to flow across magnetic field lines. These 
electric fields stem from, for example, v x B terms in Ohm's law and 
represent the electric field in a frame that is at rest with respect to 
Che plasma which is moving at velocity v . In addition, the electrons 
and ions which comprise this moving plasma respond differently when 
under the influence of magnetic and electric fields. The net result is 
that local charge separation can occur in the frame of the moving 
plasma, which, upon neutralizing, results in the generation of elec- 
trical currents. Thesa currents Chen modify Che ambient magnetic field. 
Thus, 1 simple process results in the conversion of mechanical energy 
associated with velocity fields into electrodynamic energy associated 
with the production of electric and magnetic fields. These processes 
.ir- i t: Uuli’.rly important in the photosphere and below where a large 
reservoir of mechanical energy exists (cf. III). Therefore, to under- 
stand the generation of electrodynamic energy it is essential that ob- 
servations focus upon both the power spectrum and polarization of 
velocity fields With respect to the magnetic field. 

IV. 1.2. propagation of electrodynamic energy Under certain conditions 
to be discussed below, electrodynamic energy can propagate away from 

pocceD'-'^G PfiGF WOT FILMED 
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ics $ita of generacion represented by an electrodynamic Poynting flux. 
The destiny of this energy also depends upon a number of conditions i 
but it is possible to isolate two possible courses of evolution. 
Specifically, these two possibilities are irreversible dissipation into 
therwodynamic endproducts and in-situ storage as magnetic energy. In 
this regard, there are also two different ways that the solar atmo- 
sphere can be coupled to the photospheric driver. On the one hand, the 
atmosphere could be driven in a time-averaged steady state — the 
response depending actively upon the details of the photospheric driver 
(e.g. the power spectrum and polarization of the photospheric velocity 
fields). On the other hand, however, the response might depend only 
passively upon the details of the driver. For example, in the case that 
the phtospheric driver could trigger a transient unloading of stored 
magnetic energy, These points are readily illustrated by Poynting 's 
theorem which is given by; 


J.E = - £-v.(ex3)— L 2 ^^ 


(15) 


where the first term on the right hand side represents the net flow of 
energy into cr out of a control volume (i.e., the divergence of the 
Poynting flux) and where the second term represents a gain or loss of 
’ jAdcic energy in this volume (Sec also equation 1 6 and follow- 
ing discussion). The first term clearly represents a coupling between 
regions of electrodynamic energy generation such as the photospheric 
driver and external region such as the corona. As mentioned above, this 
coupling could very well be in a time-averaged steady state. The second 
term on the right hand side of equation (15) is explicitly time depend- 
ent and thus represents the storage of magnetic energy within the 


corona. 
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In order chat Che corona ba driven in a tima-averaged steady state, 
it is essential that the generated electric fields be at least partially 
in phase with and have a component that is perpendicular to the gene- 
.’ivvl -tvnusic iiul... c.u»ai :ieUa> utcm from v < terms in 

Ol\m's law, a atmdy ataS# emvgu "input ia alukxya aaaoaiatad with maa 
flcwa aoi'oas mayn&tia fistd Zinea, both at the photospheric driver and 
within the solar corona. It is important to note that this is a 
fundamental law of magnetohydrodynamics and is completely model in- 
dependent. 


IV. 1.2. DISSIPATION OF ELECTRODYNAMIC ENERGY The third aspect oC 
electrodynamiu coupling of the photosphere, chromosphere and corona is 
che dissipation of electrodynamic energy. lonson et al (1980) point out 
that this is conveniently illustrated by the following equation for 
j.E , which is obtained by using Ohm's law, che one-fluid momentum 
cquasio.n and Maxwell’s aquations: 




(16) 




m;'. cti c'ac right hand side of equation (16) represents elec- 


tron "Joule" heating by cross-field currents. The generalized resist- 
ance ’^electric anomalous as well as classical transport 

processes. The second term also represents electron "Joule" heating, 
but in this case it results from the dissipation of field-aligned 
(i.e. force-free) currents. Here too, includes both 
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anomalous and classical transport phanotnona. Tho third term represents 
viscous heating and mechanical nccelaration of the ions due to both 
anomalous and classical processes (i.e. \^c}janical these 

tiit' ; , Ir, tlvt-j ;nie vIhccui I in halting results from 

sheared ion velocity fields that are driven by Lorentz forces. 


IV. 2. Heating by MagnetohydeodynaniQ "Hody'^ tiavea 

Magnetohydrodynamic body waves are well defined nornal mode 
oscillations which transport energy via the intrinsic elasticity and 
tensile properties of the magnetoplasma. These waves are character- 
istic of a uniform plasma and thus the addition of plasma non- 
uniformities plays only an extrinsic role (e.g. reflection and 
refraction phenomena) . There is also a second-class of waves which 
exist only in the presence of non-uniformities. These so-called 
"surface waves" will be discussed in section IV. 3. With regard 
to hydromagnetic body waves there are three different wave modes, 

* . fvSn, fast and slow modes, whose dispersion relations are 

given by: 


( 17 ) 



where 9 is the angle between k and B . 

In equation (17) u)^ represents the eigenfrequencjr of shear AlfvSn 

waves which are independent of magnetic and thermal compression and 
depend solely upon the tensile properties of the magnetoplasma. The 
eigenfrequencies (Up and (U^ represent "fast modes" in which the 
magnetic and thermal compressions are in phase and "slow modes" in which 
the magnetic and thermal compressions oppose one another. 

Before the advent of ATM instrumentation on SKYLAB and the OSO series 
of satellites, the solar corona and in particular coronal active regions 
observed to be hot but lacking in any internal fine structure. 

Thus, the theorist i^as free to treat an 
unstructured atmosphere. As 

a;.,', c.'vu-r investigations focused upon hydromagnetic body waves. 
Motivated by the qualitative expectations of Alfven (1947), Giovanelli 
(1949), Van de Hulst (1953), Piddington (1956) and De Jager (1959), 
Osterbrock (1961) provided the first quantitative study of the 
generation, propagation and dissipation of hydromagnetic body waves in a 
structureless, plane-parallel solar atmosphere. The essence of Oster- 
brock’s model was that fast-modes are generated in the highly turbulent 
convection zone (see section III) . These waves propagate upwards into 



th* chroraospher® wh«ire stetpening into ihocKu occurs just like the 
ncoustic waves. Pecause of thoir small scale size, chese. shock waves are 
efficiently dissipated within the chtomosphere and play an important 
role in the beating process, Some of the fast modes, however, may bleed 
into Che corona. Upon entering the corona, these fast modes couple to 
pure AlfvSn waves which, because they are non-compressive, are weakly 
damped and therefore have quite large damping lengths. The problem with 
this scenario is that it depends quite strongly upon Che details of the 
ambient plane-parallel chromosphere and corona. Specifically, the 
variation in the plasma $ with height above the convection zone 
critically determines the applicability of Osterbrock's approach which 
invokes geometric wave optics. Thus Che meric in Osterbrock's work is in 
his discussion of relevant dissipation processes. In particular, OsCer- 
brocW focused upon oolliHonal dissipation processes (i.e, > 

^collision^ for both 6 > 1 and (J < I regions. For example, the 
damping lengths for Alfven, fast and slow modes in a 3 > 1 region of 
Che atmosnhere are given by: 
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wher« X " n_/n^ the ioniiacion degra* and n • the viscosicy as 
n • vise * 

*..?.d in «»r*-’op Ct.l.. 'j is the wave angular frequency and is 

the collision time Cor neutral particles . 


For B < I they are given by 


= rr' .4-' 1 "^ 

L Joule eiMctien J ^ 


(I9a) 


Lp’-- 

^ S 

Lc. ^ , 

s 


(J9b) 

(I9c) 


where S is the shock strength and the shock wave period. 

Clearly these damning lengths depend strongly upon the model atmosphere, 
particularly the magnetic field and the wave period. Thus, we 
again stress the importance of determining the power spectrum and 
polaciaation of velocity fields with respect to the magnetic field as 
well as a determination of the plasma B as a function of height in the 
. 'r itfjiphere , both from the point of view of determining damping 
lengths and cransmission coefficients. 

As a foundation upon which one can motivate future inquiries, it is 
possible to establish a range of wave periods that are relevant 

to the heating of the chromosphere and corona. Noting from equations 
(I8) and (I9) that the dominant damping mechanism for Alfven waves is 
Joule dissipation whereas for fast and slow modes it is dissipation 
through shock formation, we find that '^„gyg must satisfy 
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T > 3 X lO'' B *«c for Alfvin wovti and, t > 10* v.“* for 

w chrom w " 

fast wavas otherwise they will be damped before entering 

the corona. In obtaining this constraint we have assumed that the 
scale size of the chromosphere ^chjon,osnhero ^ x jo* cm with a 
maximum density Pch^oroosphere * ^ * gm/cm' and electrical resist- 
ivity n . . ■ 3 X 10“** sec . Although we can only guess as 

^ chrofopsphttrft 

to the value of the chromospheric magnetic field, one can make reason- 
able estimates for both the quiet and the active sun. For example, if we 
assume that - 10 gauss in a quiet region and - 300 gauss 

in an active region, it follows that only Alfv^n waves with > 16 min 
couple to the quiet corona whereas waves with > 5.7 sec couple to an 
active region corona. 

Despite these constraints, Habbal et al (1979) have investigated the 

possibility of heating active regions by fast-modes with “ 2 sec . 

Although such waves arc probably dissipated within the chromosphere, 

Habbal et al’s approach might apply for longer period waves provided they 

assume a aotZiaxonat damping process. Habbal et al assume collisionless 

transit-time damping which is probably a poor assumption since the 

collision time is small compared to the wave periods of interest (i.e. 

t . ■ 1 sec < 5.7 sec with 5.7 seconds being the smallest wave 

collision 

cin penetrate into the corona). 

Iherc were two major considerations that motivated subsequent work on 
hydromagnetic body wave heating of the solar corona. First, there is the 
problem of invoking the 17KB approximation under conditions in which it 
is not applicable and second, there is the problem of very large damping 
lengths and thus apparently negligible heating of the coronal plasma by 
those waves that do manage to penetrate the chromospheric filter. By 
far, the problem that has received the most attention since Osterbrock's 
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(1961) initial work has btan that associatad with mtchanisms which in- 
crtasa tha damping o£ Alfvin wavas, tharaby decraasing tha coronal 
damping length. Tha moti\rating factor which led to thasa studies was 
tha discovery that the corona is in fact quite structured. Thus our 
previous concept of A corona with a single scale •sslio of the order of 
10*^ cm was completely altered and it became necessary to investigate 
processes that would result in efficient heating of significantly 
smaller regions such as low lying coronal loops. Despite the structured 
nature of the solar atmosphere, many authors neglected variations of 
the plasma parameters in a direction that is perpendicular to the 
magnetic field. As such, they assumed that Alfv<n waves retained their 
plane^wave nature (i.e. no cross-field variations in the wave proper- 
ties). This, however, is a reasonable assumption provided that the k'KB 
approximation is valid. 

The problem of studying the propagation of hydromagnetic waves with- 
out recourse to the eikonal approximatl.on was not addressed until after 
Belcher and Davis (1971) observed what appeared to be Al£v6n waves 
propagating in the solar wind with periods of the order of hours. Holl- 

weg (1972, 197S, 1980) recognized the important implications 
of such observations and was the first to develop a model of Alfven wave 
crc'vij’ uiou In the solar atmosphere without the restriction of the 17KB 
approA-t'acion. Although his analysis disregards various damping mecha- 
nisms, it does Indicate that long period Alfven waves of the order of 
hours can propagate through the chromosphere with sufficient power 
(P ■ 10** - 10® ergs cm”^ sec“^) to strongly affect the solar wind and 

heat at least the quiet solar corona. Subsequently, Wentzel (1974, 1976), 
Uchida and Kaburaki (1974) and Melrose (1977) claimed that Alfven waves 
can be converted into fast-mode waves in the corona from bends in the 
coronal magnetic field. If the scale length of these bends is of the 
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order o£ the wavelength of the ALfvCn wave then there results a ^ 
induced density compression which is characteristic of a fast mode. 
This density compression repidly leads to steepening of the wavefront 
into a shock, and subsequent dissipation in about one wavelength. In 
p.*iucipl« this seems raasonaole, but there is a fundamental problem as 
to how energy can be coupled into the electron and ion microstructure 
of the plasma. Specifically, the conversion of Alfvfin waves into 
fast modes occurs only when the wavelength of the Alfvln wave is of 
the order of the loop size. As such, the resulting shock thickness is 
also of the order of the scale size of the coronal loop and one must 
be quite careful when investigating heating of the electron and ion 
populationis. 

Wenezel (1974, 1976) also Invoked weak non*linearicies into hts 
analysis of AlfvSn waves. Although AlfvSn waves are non-compressiva to 
first order in the wave amplitude, they are to second order com- 
pressive. Specifically, the second order pressure fluctuation 
Pj * 4ir^ / 8 it (Bj is the first order magnetic fluctuation)” 
and the second order velocity displacement v^ ■ 4 tt^ v^ / B^^ 

« ambient field strength) results in a propagating density wave 
whose intensity is given by P^v^ . Assuming that this density com- 
. ?ouples energy into heat, after travelling about one wave- 

A » the damping length associated with the driving 

AlfvSn wave is 

‘■damping ■ (B./P.) (v, . 

It should be noted that 


cat 



35 . 


v;«ntz«l coniidartcl only thoi« cas<it in which tht WKB approach was a 
good sppcoxiraation (i.«. waveUngth small comparsd to tht magnttic 
fltld't radius of curvatuca). 

From Wtncztl's analysis it it possiblt to dtttmint what wavt 
^a vs.b? that his r.-icaaniam rasult In o£?icient 

hasting ot tha corona. Noting that noting that 

tht damping langth roust ba tmallar than tha charactaristic seals size 

of tha corona I ^ - I x lo“ cm , it immadiataly follows from his 

corona 

axprassion for tha damping length 

must satisfy < 3.9 x lo‘* ^ 

it tha maximum observed velocity ditplaeament in the solar corona . 

Since Vj®*'* » 2 x |0* cm/sac (Barring and Feldman, 1974 ), we find 
chat only waves with x^ < 24 min can efficiently heat quiet regions 
whereat only waves with X^ < 5.7 sec can heat active regions. Although 
Wentzel's damping mechanism is consistent with the heating of quiet 
regions there is some question at to its applicability in active regions 
since only waves with X^ > 5.7 sec escape from chromospheric damping. 
It should also be pointed out that despite the large wavelengths 
asj 50 cl*;cud wicn wa>/« periods < 24 rain (i.e. X < 3 x jo^® cm) 
Wentzel's analysis is still approximately valid since quiet regions are 
-.r.h much less structured than active regions. It would, however, 

■s~. ■' iCavi'C'S to davalop a model of quiet region heating that in- 

corporates the virtues of both Hollweg's 0978) non-eikonal method of 
studying long period AlfvSn waves and Wentzel's (1974, 1976) non linear 
damping mechanism. 


3. Heating by Magnetohydrodynamia ''Surfaae*' Waves 


lonson (1977, 1978, 1980) was the first to recognize the importance 


PAGE IS 

"mi QVAun 



36 . 


of global afruccura oscillations in chs inhscsntly structured active 
region corona, lonson, following Chen and Koscgawa (19/4) and Hasegawa 
and Chan (I976)» has stressed that the response of a plasma structure 
to external driving is qmliiativi*ly difftmnt from that which occurs 
• 4 . ...t a..a» ru'.wurui swvtLug. this difference jt^vsi from the fact that a 
plasma structure hosts a continuous spectrum of internal oscillators 
(i.e. hydromagnetic body waves), each representing the oscillatory 
characteristics of an inf initessimally small piece of the plasma. For 
example, consider a coronal loop which is characterized by a parabolic 
depression in the Alfvtn speed across the minor radius of the loop 
(i.e. v^ ■ * given driving period tp , there will be a 

continuous spectrum of excited Alfvin waves whose wavelengths satisfy 
X(x) ■ ^ existence of such a continuous spectrum of Alfvin 

waves (other body waves could be excited as well) results in both a 
''global" and "local" response of a nonuniform plasma structure to ex- 
ternal driving. The global response reflects the integrated response of 
all the internal oscillators which in this case are Alfvin waves. As 
such, it is associated with a spatially coherent oscillation of the 
entir-* plasma »cruc£ure. fhe local response, on the other hand, 
corresponds to that of a single internal Alfvin wave and is essentially 
" ' ‘»!4 fiuld line oscillation, lonson (1977, 1978) originally called 

-V.4* i. civil ■structure oscillations Alfvinic mir'faae waves to stress 
chat the presence of noymn'ifovmity is essential to their existence 
(c.f. Wentzel 1978, who discusses other surface waves). It is im- 
portant to note that a discontinuous surface is not required. The wave- 
length of the excited surface wave is given by 



* 
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and in os««nc« raprtitnCf th« Cicld-alignid auco-correlacion l«ngch 
aisociattd with Cht flow of cltccrodynamic tn«rgy. lonson (1977, 1978) 
has also pointed out that: standing surfnca WAvas could rtsult providad 
■ n%it whiro K is the loop length and n is some integer. In 
addition, Zonson (1980) has pointed out that other global structure 
oscillations such as raagnetosonic cavity modes exist in the natural 
cavities associated with regions between magnetosonic wave reflection 
points. In particular, multiply reflected magnetosonic body waves phase 
mix and establish a discrete spectrum of standing (in the direction of 
the non-uniformity) magnetosonic cavity modes. In this regard, Zonson 
(1980) has proposed that coronal structures be interpreted as either 
hydifomagnetia vesonanaa oavibiea such as coronal loops or hydi*magn3tia 
mv9gu.id«8 such as coronal holes. This viewpoint has also been dis- 
cussed by Qucrfcld and Hollweg (1980). 

One of the most interesting aspects of global structure oscillations 
such as Alfvinic surface waves and magnetosonic cavity modes is their 
linear coupling to the Alfvgn wave continuum at local Alfvdn resonances 

where for example A(x*x ^ ) ■ X. . This feature was noted by 

^ resonance A 

lonson (1977, 1978) in connection with the linear degradation of Alf- 
vSnic surface wave energy into localised body wave energy at a spatial 
*... .i vitHln n coronal loop. The so-called process of "resonance 
absorption" (lonson 1978, and references therein) is a specific, but 
extremely important example that has been found to control the transfer 
of energy from an external MHD "driver" to the interior of the "driven" 
plasma structure. In simple, qualitative terms, this process is 
associated with the existence of an Alfvenic resonance within the 
plasma structure where the local impedence (i.e. related to the local 
Alfven frequency) matches chat of the external driver (i.e, related- to 
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thtt oicillacion frequency of the spatially coherent external driver). 

At the spatial resonance, where impedance matching occurs, dissipation* 
processes, such as ion viscosity, play a critical role in converting 
the incoming Poynting flux of electrodynaraic energy into heating the 
co-’caal plaotiui. Since Lhu dissipation race iaereatisiS aich decreasing 
scale size and since the resonant absorption layer is naturally small 
compared to the global system, it follows that plasma heating and 
acceleration is localized within Che resonance region. In this regard, 

Che resonance absorption process facilitates the dissipation of large 
scale mhd disturbances by providing a "sink" (i.e. the resonance 
absorption layer) for large scale electrodynomic energy which, because 
of its inherently small spatial scale, results in an efficient con- 
version of this energy into thermodynamic endproducts. 

It is important to note that only within the resonance absorption 
layer, will the effects of local irreversibilities make an explicit 
showing. That is to say, they will define the spatial extent of Che 
resonance absorption region and determine the local amplitude of the 
fluctuating electrodynamic fields. 

Ihus, the relevance or lonson's work is in his pointing out the exist- 
ence of large scale structure oscillations such as the AlfvSnic 

-rav.is and his discussion of how the energy associated with 
thw-j’ larsc scale oscillations is coupled into the microstructure of 
the plasma, i.e., the resonance absorption process. Although lonson did 
not self-consistently include all aspects of generation, propagation 
and dissipation of electrodynamic energy into his analysis, he did 
underline the potential importance of 300 second oscillations in heating 
solar coronal loops. 
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IV. 4. Mating by Cut'Mnia 

IV, 4.1. CORONAL HEATING BY ANOmOUS JOULE DISSIPATION Over Cha Xait 
faw yaari numerous proposals have been made suggesting that coronal 
U*ieL'U v.'S -.lt from she dissipation of currents (Tuc'ter 1971; Nolte 
et al 1977; Rosner at al I97S; Hlnata 1979, 1980). There are basically 
three physical mechanisms capable of dissipating such currents: Joule 
dissipation, reconnection and double layers. Joule dissipation and 
reconnection are closely related, differing only in that pure Joule 
heating causes no topological changes in the magnetic flUK surfaces of 
Che magnetic structure, vhile reconnection causes such topological 
changes together with Joule heating and strong convective flows. Here 
we will review critically the possibility that coronal heating tsay 
result from either Joule heating or reconnection since double layers 
have not, as yet, been proposed as a coronal heating mechanism (for a 
review of these three mechanisms c.f. Spicer and Brown 1980). 

To develop a theoretical framework in which we can study Che 
feasibility of either Joule heating or reconnection as viable coronal 
heating mechanisms we must specify the physical constraints on such 
heating mechanisms imposed by both theory and observations. The two 
most obvious constraints are chat in steady statai (i) the local 
- l.ri input from the heating mechanism ^ (the local 
losses) and (ii) the power input from the heating mechanism 
integrated over the fractional volume of the magnetic structure that 
is heated must roughly equal the integrated power loss from the entire 
volume of Che magnetic structure. 

aV V 




( 21 ) 
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where AV is the fractional volume in which the heating mechanism 
operates and V is the total volume of the structure. Obviously if 
heating is occurring uniformly throughout the volume and at 

enrh poim. then AV ■ V . However, both reconnection and Joule heating 
ru.iuire ch.u: ».V/V «.s i so that ‘ these re- 

quirements fundamentally alter the energy balance arguments previously 
proposed (e.g. Rosner et al 1978b). 

To demonstrate that both Joule heating and reconnection require 
AV/V « 1 and thus >> we note that at the very laast the 
heating mechanism must overwhelm radiation losses, ■ n* f(T) (see 
section II). 

For irreversible heating both Joule heating and reconnection result 

in £„ • n J^ , where the resistivity n “ V . O) , v . being the 
H ’ ' ei pe ' ei 

effective electron-ion collision frequency, the plasma frequency, 

J ■ - ne ^ and being the drift velocity of the electrons relative 

to Che ions. Using £j^ ~ we find 

Vjf > (22) 

Taking ( JJ) as an equality and substituting both the expressions for Ej^ 

and £„ into ( 22) we find / n^ f(T) d V « / n^ f(T) d V which 

“ AV V 

■ ^ -satisfied if AV •« V . However, by using Ampere's equation 

' A v-'-/o) J we find that the characteristic width of the current 

channel satisfying (22) is 


S£ 


< 


c 

ht.2 
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where 6B is the change in the current produced magnetic field across 
the current channel. If we take the heating to occur along the whole 
length L of a cylinder of radius r^ then the ratio of AV/V ^^/r^ , 
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which is usually isuch lass than sue. To make AV/V * I would cause the 
magnetic field associated with the current density required by (22) 
be ~ 10^ g in the corona. Hence, we arrive at a contradiction so that 
''•s I and thus >>6^ from which we must conclude that 
>> f(T)/ne" . However, there is a limit to how high can become 
before exciting various micro instabilities that are capable of in- 
creasi'ig n by several orders of magnitude, t^en reaches a 
magnitude roughly comparable to or greater than the ion sound speed 
C ■ (k- T /m.)^ these instabilities are excited and one has to take 

S o € I 

into account the phenomenon called "anomalous resistivity" (section 
II. 3). Under these circumstances the Joule heating becomes anomalously 
large so that >> is easily satisfied. Anomalous Joule heating 
is the basis for the coronal heating models proposed by Tucker (1973) 
and Rosner et al (1978) and examined in greater detail by Hinata (1979, 
1980) and Vlahos (1979). 

We stress that an ideal force free field is incapable of supporting 
such an anomalous heating since for more than a dissipation time force 
free fields have zero Poynting vectors and thus (Stt)”* 3B^/3t ■ - 
(see eq. (1)) . Using £ “ (Ohm's Law) and £ ■ (ac/Air) ^ for a 

force free field, it follows that ■ B^^ exp (-t/t^j) where 

* ^ and a ■ V X B/£ . Thus if the current associated with 
a force-free field is to drive unstable the ion-acoustic instability 
(or any other instability that is capable of causing anomalous resist- 
ivity) it must have an a ** \/6Z as defined by ( 23) on some field line 
so that the current will only last for a short time. Hence, coronal 
heating due to anomalous Joule dissipation and driven by a force-free 
field current is a one shot heating process so that steady state heat- 
ing by force-free field currents is impossible. Thus, for any semblance 
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of sceady stnCe heating to be possible we must have a non'-zcro Poynting 
vector, Furtheri we can rule out a Jj| associated with an equilibrium 

, since in equilibrium is determined by the pressure gradient 
which in turn would necessarily be extremely excessive for coronal 
condi cions if ic to support a sufficient to drive a current un- 
stable to microscopic instabilities. Hence, if is to drive the 
current unstable to any instability that causes anomalous resistivity, 
convective flows must exist in the coronal structures that are 
generating . 

Let us assume that such flows do in fact exist so that some 
semblance of steady state anomalous Joule heating can be achieved. 
Physically this implies that the drift velocity of the current must 
satisfy * Cg and the instability will be near marginal stability. 
Marginal stability is simply a statement of the fact that microscopic 
mechanisms such as the ion acoustic instability grow very rapidly. The 
ion acoustic instability growth time is tg ” ^pe”* 

saturates in ^ 20 tg so that if one insists that. such mechanisms 
play important roles in coronal heating the instability roust be con- 
stantly driven by external means towards instability. Thus for a 
quasi-steady state to exist the threshold conditions must be marginally 
■ ‘.•T.iheimar and Boris 1977; Spicer and Brown 1980). A 

□argiaaiiy stable system can pulsate about some quasi-steady state, if 
the dissipation time is shorter than the build up time. The heat 
pulses have a temporal width given by the dissipation rime while tne 
pulses will be temporally separated hy che build up time. However, if 
the build up time is shorter than the heating time the heating will 
continue until the dissipation time equals the current hm'ld up time 
because the threshold condition Cg will increase so Chat will 
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decrease until the dissipation time equals the build up time. At this 
point a quasi-steady state will be achieved with the instability rapidly 
switching on and off. Averaged over the seeing time of the average solar 
instrumant such a coronal heating mechanism would appear in steady 
state, so that the assumption that the instability is at marginal 
stability is a good one if such a mechanism is in fact operating to heat 
coronal structures. 


IV. 4, 2. HE.VTING BY RECONNECTION Heating of the corona by re- 
connection has been proposed by Levine (1973), where he argued that if 
large nund).ers of ^ driven neutral sheets were to exist in the solar 
atmosphere the Joule heating and particles accelerated by the reconnect- 
ion process would serve to heat the corona. His arguments were 
qualitative in nature, with little support either observationally or 
from detailed theoretical calculation. Since Levine's hypothesis is 
basically an application of a flare mechanism to coronal heating, and 
Rosner ct al (I97Sb) have resurrected the AlfvGn-Carlqvist flare Modal 
(Alfven and Garlqvist 1967; Smith and Priest 1972) as a coronal heating 
•rode'. ’! t would be useful to examine briefly those flare models that use 
r.i :o»nec tion by tearing modes (Spicer 1976, 1977; Colgate 
1978) to see if they would be practical as scaled down coronal heating 
models . 

Reconnection either by neutral sheets (c.f. e.g. Priest 1976) or by 
tearing instabilities (Furth, Killeen and Rosenbluth 1963; hereafter 
FKR) represents a mechanism by which energy stored globally is released 
locally in a concentrated form. Depending on the collisionality of the 
plasma and the characteristic magnetic field gradients the reconnection 
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process can be collisional » oollisionless '^gi) 

or send-collisional (Yj^j * ) , where and v^. are respect- 

ively Che tearing mode growth rate and the electron-ion collision 
frequency. In practice if Che plasma is collisionless or semi- 
collisional the dicficulcies associated with anomalous Joulu healing 
again arise: that is, extremely steep field gradients and a strong 
e:<ternal driver maintaining Jy are again required. On Che ocher hand, 
the constraints on collisional tearing are not as severe. Nevertheless, 
regardless of what the collisionaliCy of Che reconnection process is 
there are a number of general arguments that can give us some insight 
into the viability of reconnection by tearing modes as a coronal heat- 
ing mechanism. 

Reconnection in neutral sheets is treated in steady state, while in 
more complex magnetic configurations such as loops the tearing mode, a 
dynamic mechanism, is responsible for the reconnection process. The 
tearing mode gives rise Co inductive electric fields parallel to the 
magnetic field, and it is this electric field that is reponsible for 
either Joule heating and weak particle acceleration (collisional tear- 
ing) or strong particle acceleration (semi-collisional and collision- 
less tearing). This parallel electric field, 6E|| , is generated within 
i whose thickness A is small compared to the characteristic 

scale sices or the magnetic field configuration . It is within 
this layer that most of the dynamics of the instability takes place. 

Tne first point to remember about the tearing instability is that its 
wavelength parallel to the magnetic field is infinite, i.e. 

A|j ■ * (k . ^ ■ 0) , while its wavelength perpendicular to ^ must 
satisfy 6^/Aj^ < I (FKR 1963) . The fact that Xjj ■ “ implies that a 
significant spatial scale of the magnetic configuration parallel to B 
is involved in the reconnection process, while a spatial scale 
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psrpendicular to Z with a seals sise t 51 will be involved. 
Hence, the tearing mode will manifest itself as a long filamentary 
brightening parallel to ^ with a thickness ^ Xj^ . As a general rule 
the smaller the ratio Sl/X^ the more effective the tearing mode will 
be in extracting magnetic energy t'roci global storage and releasing it 
within A (A « 6S) . This follows because wavelengths that satisfy 
5-t/Xj^ S I strongly bend the magnetic field and thereby generate the 
largest restoring forces. These large restoring forces in turn lead to 
very low amplitudes for the modes at saturation. Thus little energy 
will be extracted from the magnetic field. On the other hand those 
modes that satisfy 6t/Xj^ « 1 generate the weakest restoring forces 
and grow to large amplitudes before saturation, thereby extracting a 
larger fraction of the available global free energy supply. However, 
long wavelength modes not only lead to significant Joule heating (or 
particle acceleration) but also to strong convection and therefore 
strong disruptions of the magnetic configuration. Our arguments imply 
that short wavelength modes (5t/Xj^ ~ I) will not significantly heat 
the plasma while the long wavelength modes (fS/Xj^ « I) can heat the 
plasma but will also cause large scale disruptions of the plasma- 
magnetic field configuration. 

ire two types of tearing modes treated by FKR (1963), the 
alcv-Ccacins modes and the fast-tearing nodes. As a rule the thinner 
A relative to Xj^ the more likely that the mode will be a fast tear- 
ing mode. The non-linear slow-tearing modes have been treated by 
Rutherford (1973) and were shown to saturate at low amplitudes, thereby 
extracting little magnetic energy. However, the fast tearing mode is 
known to grow exponentially right up until saturation and for this 
reason was proposed as an excellent candidate for a flare mechanism 
(Spicer 1976, 1977). More recently these modes were shown analytically 
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to bft capablft of heating the plasma (Spicer I960). However it was found 
that only after many e-foldings significant heating begins to occur 
because the perturbed parallel inductive electric field is very weak 
initially and must grow over many e-folding times before it is compar- 
'.cia t; -Uc«:ci. £Ulu driving tha jquilibrium Jj. . Since the 
slow-tearing mode saturates at low amplitudes it is highly questionable 
whether it should be considered as a good candidate for coronal heating, 
while the fast-tearing mode which does lead to heating causes strong 
convective disruptions. In cylindrical geometry these disruptions 
should appear as a kinking of the cylinder (or loop) (resistive kinks). 

Unless such behavior is occurring in quiescent coronal loops it is 
highly questionable whether tearing modes can be considered to play a 
role in the coronal heating mechanism. 

We should also add before closing this section that collisionless 
and semi collisionless tearing modes are not likely to give rise to 
scaling laws similar to those obtained by Rosner et al (1973a), because 
these scaling laws result essentially from balancing thermal collisional con- 
duction against radiation losses . However, since semi-collisionless and 
collidionless tearing modes result in supra-Cheroal electrons the 
typical mean free path of these electrons can ''xsily exceed the scale 
?is® o*' chs magnetic configuration itself invalidating the use of 

conJuction as a transport mechanism under such circumstances. 

Hence, Che resulting scaling laws are likely to be quite different. 

IV. 5. k "Unified'^ Approach to Solar Coronal Heating 

lonson et al (1980) have developed an approach 
to electrodynamic heating of the solar corona that unifies all previous- 
ly proposed mechanisms and at the same time self-consistently 
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incorporactf ftldccrodynamic tnargy gcnaracion, propagation and 
difftipation, Th«ir so-call«d '*eransml»sion lino approach” is basad upon 
traating solar plasma strucfcuras as eithar ”closad" or "open” trans- 
mission linos (o.g. coronal loops and coronal holas) charnctarizad by 
impadancasy Q-valu#s, L/R , RC and LC timas, This approach is valuable 
in that it strassas the fundamental aspects of alecnrodynamics • viz,, 
inductance, capacitance and "resistance” while at the same time 
providing a formalism which describes the self-consistent electro- 
dynamic coupling of the photosphere, chromosphere and corona. This 
coupling is described by a simple circuit equation which includes for 
currents flowing along and across the magnetic field, i.e. j|| and 
an inductance L , a generalized "resistance” R and the plasma 

t 

capacitance C « Zonson et al (1980) point out that in certain limits, 
these equations reduce to those associated with more traditional 
heating mechanisms. For example, the "quality” or Q-value of the 
circuit is defined as Q ■ (L/CR^)^ and is a measure of how good an 
oscillator the circuit is. If a particular coronal structure is 
characterized by Q » I then this structure is a good oscillator and 
the circuit equation reduces to that associated with hydromagnetic 
wave-heating theories. In this case, power absorption is highly peaked 
• r.'Eur.'il oscillation frequency of, for example, a coronal 

loop. For Che most part, the natural frequency of oscillation for most 
coronal loops appears to be centered around 300 seconds. Therefore, 
lonson et al (1980) claim that high-Q coronal loops should be treated 
as hydromagnetic waveguides. If, however, a particular coronal 
structure is characterized by Q « i , then the circuit equation 
reduces to that associated with "current-heating" theories. In this 
case, power absorption is not highly peaked about some resonance 
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£r<qutncy and indued, resonancca do not even exiit. 

lonson ct al (1980) alio point out that the above theoty ii subject 
to potentially obasi*t>abt« boundary condition! (v and B) . At such, this 
npp coach could be very important in allowing an efficient ai*imilation 
of data from ongoing and future solar missions such as the Solar 
Ma)(imua Mission and Spacelab missions. 
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V. CONCLUSION 

Iht mtchanism oi the heacing of tho lolac chromosphora and corona 
remaina an opan quaaeion aapacially if it coaiaa to axplatn datailad 
a:n9»»r»h?ri4 j^cucturas. Th3 «n;>rtnln*:y i» aaaontiaUy ilua to our 
poor knowladga of tha magnacic fiald and tha velocity fiald apactrum 
with raapact to tha magnacic fiald. Xha magnatic fiald in tha chrono- 
aphara and corona appaara to ba known only qualitativaly, whila for a 
haating theory a rathar datailad quantitative knowledge of the 
magnetic fiald ia required. Wa have amphaaizad chat it ia mainly the 
value of 0 whoaa diatribution throughout tha aCmoaphare ahould ba 
known. According to whathar 0 » I or 0 « I wa have divided tha 
existing thaoriea in 

a) the acoustical thaorias, such as tha well<-known shock wave theory 
which has bean widaly applied also to explain other stellar 
eoronae; 

b) the alactrodynamic theories , .which may ba divided into heating by 
body HHD waves , surface MKD waves and currant haating. 

As a conclusion it can be stvited tint Che acoustic theory is losing 
its range of applicability concerning the expected distribution of 0 . 
Figure I schematically indicates the range of validity of the acoustic 
. ... elect rodynamic theories. 

The acoustic theory might be correct with little modifications for 
most of the solar chromosphere but as soon as Che atmosphere shows a 
high degree of structure such as in the corona and transition layer 
Che magnetic field must play a dominant role. 

As far as the electrodynamic theories are concerned it appears that 
also the current heating theories have a small range of applicability. 
They are particularly important for flare heating and for flare - like 
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coronal phcnofflcna. 

Ho«c promliing are eha wave eheorita. Bodjr &$ well at surface waves 
appear Co play an iBiportant role. In the "open” magnetic structures 
body waves may be crapped and these structures may ace as a kind of 
• ;v- i'l* v.vvoa. The ■na-jn.jcic 3rructu**cs may act as 

resonant cavities for Alfvfnic surface waves and may therefore act as 
a special kind of absorber for Alfvin waves. We have stressed the 
importance for a unifying and consistent theory of the formation and 
heating of coronal structures that eventually bridges the gap between 
Che photospheric "driver" (e.g. the velocity and magnetic fields as 
they should be observed in great detail in the photosphere) and the 

observed chromospheric and coronal radiation structures. A first 

* » 

attempt in this direction in which coronal structures are considered 
as electrodynamic circuits seems to le.?^ to a promising 
and original way of describing the solar outer atmosphere. 

Finally the reader is warned not to reject the acoustic theory 
categorically, because it is "en vogue" to do so presently, but to 
consider Che ultimate electrodynamic theory as a more powerful theory 
that weutually includes the acoustic aspects of chromospheric and 
coronal heating. The prime question which remains unanswered is 
namely whether the corona is heated by essentially longitudinal waves 
. - ■'r.Kcic ’wves and magneto-acoustic waves (fast and slow 

modes) or by essentially transversal waves such as the AlfvSn waves 
(body and surface waves). 

Due to the paucity of the observations of v and B in the solar 
atmosphere this question can not be settled. The only facts that seem 
to be generally accepted is that the lower atmospheric oscillations 
are acoustic in nature and the magnetic fields are strong and con- 
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c«ner«t«d in flux tubsi in the photoiphers. Until we hnv« observed the 
etate of motion of the chromoiphere and corona with reepect to the 
maanetic fiulw we have to live with two type* of theory, tho acoustic 
theory and the more general electrodynomic theory. 
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Caption to tht figura. 


60 


Schtroatic rapresentation of the magnetic atructura of the solar 
atmosphere. The scale is highly arbitrary. The lower atmosphere 
where to a large extent $ « I is mechanically coupled to the 
I’or.vective nnd turbulent motions and the 300 sec. oscillations. 
The corona where g « I is electrodynamically coupled to the 
lower atmospheric velocity field. In the suhcoronal layers it 
is the mechanical energy flux and in the corona the Poynting 
flux that is responsible for the energy transfer. 



